------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\Stijn\Documents\_Projects_\Research\Papers\Prestatieanalyse\Persistence_in_Top\Analyse\130401_ReStat_scripts_data\su
> bmitted_to_ReStat\citanalysis.log
  log type:  text
 opened on:   8 Apr 2013, 13:23:19

. 
. use citdata.dta, clear

. 
. stset end, origin(time max(zapdate,d(01oct1992))) failure(h2y==1) id(persnr) time0(begin) scale(365.25) exit(finalobs==1)

                id:  persnr
     failure event:  h2y == 1
obs. time interval:  (begin, end]
 exit on or before:  finalobs==1
    t for analysis:  (time-origin)/365.25
            origin:  time max(zapdate,d(01oct1992))

------------------------------------------------------------------------------
     7764  total obs.
     1034  obs. begin on or after exit
------------------------------------------------------------------------------
     6730  obs. remaining, representing
     1034  subjects
      776  failures in multiple failure-per-subject data
 6729.604  total analysis time at risk, at risk from t =         0
                             earliest observed entry t =         0
                                  last observed exit t =  8.999316

. replace _t0=round(_t0)
(4483 real changes made)

. replace _t=round(_t)
(5375 real changes made)

. 
. stcox male age maind_B maind_C maind_E maind_G maind_H maind_I maind_M maind_N maind_P maind_R maind_Z fac_wet fac_twt fac_lbw fac_gen f
> ac_far r_doc1 r_hdoc1 r_hl1 r_rest1 yrsrank headuni1 zapaft92 yrszap tchload cumGOAlag1 cumOTlag1 avgco avgcom fulltiun if prefirst2y==1
>  & pubm_pre92h==0, robust cluster(persnr) /* first top: table 3, dependent variable = citations */

         failure _d:  h2y == 1
   analysis time _t:  (end-origin)/365.25
             origin:  time max(zapdate,d(01oct1992))
  exit on or before:  finalobs==1
                 id:  persnr

Iteration 0:   log pseudolikelihood =  -681.9746
Iteration 1:   log pseudolikelihood = -631.86572
Iteration 2:   log pseudolikelihood = -625.57473
Iteration 3:   log pseudolikelihood = -625.23536
Iteration 4:   log pseudolikelihood = -625.19458
Iteration 5:   log pseudolikelihood = -625.18019
Iteration 6:   log pseudolikelihood = -625.17489
Iteration 7:   log pseudolikelihood = -625.17295
Iteration 8:   log pseudolikelihood = -625.17223
Iteration 9:   log pseudolikelihood = -625.17197
Iteration 10:  log pseudolikelihood = -625.17187
Iteration 11:  log pseudolikelihood = -625.17183
Iteration 12:  log pseudolikelihood = -625.17182
Iteration 13:  log pseudolikelihood = -625.17182
Iteration 14:  log pseudolikelihood = -625.17181
Iteration 15:  log pseudolikelihood = -625.17181
Iteration 16:  log pseudolikelihood = -625.17181
Iteration 17:  log pseudolikelihood = -625.17181
Iteration 18:  log pseudolikelihood = -625.17181
Iteration 19:  log pseudolikelihood = -625.17181
Iteration 20:  log pseudolikelihood = -625.17181
Iteration 21:  log pseudolikelihood = -625.17181
Iteration 22:  log pseudolikelihood = -625.17181
Iteration 23:  log pseudolikelihood = -625.17181
Iteration 24:  log pseudolikelihood = -625.17181
Iteration 25:  log pseudolikelihood = -625.17181
Iteration 26:  log pseudolikelihood = -625.17181
Iteration 27:  log pseudolikelihood = -625.17181
Iteration 28:  log pseudolikelihood = -625.17181
Iteration 29:  log pseudolikelihood = -625.17181
Iteration 30:  log pseudolikelihood = -625.17181
Iteration 31:  log pseudolikelihood = -625.17181
Iteration 32:  log pseudolikelihood = -625.17181
Iteration 33:  log pseudolikelihood = -625.17181
Refining estimates:
Iteration 0:   log pseudolikelihood = -625.17181
Iteration 1:   log pseudolikelihood = -625.17181
Iteration 2:   log pseudolikelihood = -625.17181

Cox regression -- Breslow method for ties

No. of subjects      =          806                Number of obs   =      4219
No. of failures      =          115
Time at risk         =         4219
                                                   Wald chi2(32)   =   7059.20
Log pseudolikelihood =   -625.17181                Prob > chi2     =    0.0000

                               (Std. Err. adjusted for 806 clusters in persnr)
------------------------------------------------------------------------------
             |               Robust
          _t | Haz. Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        male |   9.618537   5.693352     3.82   0.000     3.014883    30.68652
         age |   .9198755   .0191524    -4.01   0.000     .8830931      .95819
     maind_B |   1.949488   .9107572     1.43   0.153     .7802938    4.870605
     maind_C |   1.782398   .7694439     1.34   0.181     .7648025    4.153939
     maind_E |   1.667519   .7869321     1.08   0.279      .661262     4.20502
     maind_G |   2.17e-15   1.28e-15   -57.43   0.000     6.86e-16    6.88e-15
     maind_H |   4.937699   2.404132     3.28   0.001      1.90143    12.82239
     maind_I |    1.57331   .6551898     1.09   0.276     .6955683     3.55868
     maind_M |   1.935992    .755026     1.69   0.090     .9014464    4.157836
     maind_N |   6.944715   4.248755     3.17   0.002     2.093592    23.03652
     maind_P |   2.922788   1.723155     1.82   0.069     .9203732    9.281764
     maind_R |   .8061367   .8985516    -0.19   0.847     .0907039     7.16459
     maind_Z |   .6321551   .4183459    -0.69   0.488     .1727897    2.312754
     fac_wet |   1.550161   1.289431     0.53   0.598      .303629     7.91426
     fac_twt |   3.056133   2.230431     1.53   0.126     .7310501    12.77607
     fac_lbw |   9.900769   7.256656     3.13   0.002     2.353913    41.64352
     fac_gen |    6.81003   4.363491     2.99   0.003     1.939752    23.90848
     fac_far |   1.015305   1.270158     0.01   0.990     .0874449     11.7885
      r_doc1 |   .3064223   .1205673    -3.01   0.003     .1417098    .6625838
     r_hdoc1 |    .392393   .1200427    -3.06   0.002     .2154358    .7147016
       r_hl1 |    .402822   .1479971    -2.47   0.013     .1960564    .8276473
     r_rest1 |   .3938215   .2007553    -1.83   0.068      .145007    1.069571
     yrsrank |   1.001383   .0322635     0.04   0.966     .9401032    1.066658
    headuni1 |   1.778926    .450334     2.28   0.023     1.083121    2.921723
    zapaft92 |   1.178135   .3600696     0.54   0.592     .6472102    2.144594
      yrszap |   1.006162   .0370305     0.17   0.867     .9361395    1.081422
     tchload |   .9591198   .0286096    -1.40   0.162     .9046537    1.016865
  cumGOAlag1 |   1.523335   .3291946     1.95   0.051     .9973552    2.326704
   cumOTlag1 |    1.25632   .2110068     1.36   0.174     .9039317    1.746083
       avgco |   1.242073   .0649361     4.15   0.000     1.121105    1.376095
      avgcom |   .8014097   .0419275    -4.23   0.000     .7233061    .8879472
    fulltiun |   1.818356   .7944067     1.37   0.171     .7723328     4.28108
------------------------------------------------------------------------------

.         
. stcox male age maind_B maind_C maind_E maind_G maind_H maind_I maind_M maind_N maind_P maind_R maind_Z fac_wet fac_twt fac_lbw fac_gen f
> ac_far r_doc1 r_hdoc1 r_hl1 r_rest1 yrsrank headuni1 zapaft92 yrszap tchload cumGOAlag1 cumOTlag1 avgco avgcom pubm_pre92h pasttop2y pas
> ttop2ym fulltiun, robust cluster(persnr) /* repeated top: table 3, dependent variable = citations */    

         failure _d:  h2y == 1
   analysis time _t:  (end-origin)/365.25
             origin:  time max(zapdate,d(01oct1992))
  exit on or before:  finalobs==1
                 id:  persnr

Iteration 0:   log pseudolikelihood = -3966.6936
Iteration 1:   log pseudolikelihood = -3756.9258
Iteration 2:   log pseudolikelihood = -3457.5699
Iteration 3:   log pseudolikelihood = -3418.1737
Iteration 4:   log pseudolikelihood = -3417.5983
Iteration 5:   log pseudolikelihood = -3417.5972
Refining estimates:
Iteration 0:   log pseudolikelihood = -3417.5972

Cox regression -- Breslow method for ties

No. of subjects      =          908                Number of obs   =      5284
No. of failures      =          624
Time at risk         =         5284
                                                   Wald chi2(35)   =    866.72
Log pseudolikelihood =   -3417.5972                Prob > chi2     =    0.0000

                               (Std. Err. adjusted for 908 clusters in persnr)
------------------------------------------------------------------------------
             |               Robust
          _t | Haz. Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        male |   6.461687   2.958298     4.08   0.000     2.634194    15.85054
         age |   .9459195   .0104544    -5.03   0.000     .9256496    .9666332
     maind_B |   1.539426   .4001928     1.66   0.097     .9248652    2.562354
     maind_C |   2.206036    .479627     3.64   0.000     1.440615    3.378137
     maind_E |   2.230271   .6927752     2.58   0.010     1.213261    4.099785
     maind_G |   1.555926   .4569376     1.51   0.132     .8750074    2.766726
     maind_H |   2.856211    .762606     3.93   0.000     1.692467    4.820149
     maind_I |   1.367361    .341058     1.25   0.210     .8386284    2.229444
     maind_M |   1.677309   .4042593     2.15   0.032     1.045826     2.69009
     maind_N |   1.605051   .6590713     1.15   0.249     .7177374    3.589318
     maind_P |   2.148805   .6319351     2.60   0.009     1.207452    3.824055
     maind_R |   1.970455   .6822109     1.96   0.050     .9996858    3.883913
     maind_Z |   1.314817   .4240048     0.85   0.396     .6988247    2.473788
     fac_wet |   1.608114   .9967753     0.77   0.443     .4772073    5.419093
     fac_twt |   1.474272   .9137231     0.63   0.531     .4375429    4.967463
     fac_lbw |   3.234183   1.966279     1.93   0.054      .982339    10.64799
     fac_gen |   3.048393   1.803219     1.88   0.060     .9562223    9.718136
     fac_far |   1.070557   .7058264     0.10   0.918     .2940399    3.897745
      r_doc1 |   .2684512   .0607241    -5.81   0.000     .1723147    .4182236
     r_hdoc1 |   .5340166   .0965025    -3.47   0.001      .374743    .7609847
       r_hl1 |   .6450512   .1042866    -2.71   0.007     .4698713    .8855425
     r_rest1 |   .7553383   .1598293    -1.33   0.185     .4989161    1.143551
     yrsrank |   .9969018   .0142456    -0.22   0.828     .9693683    1.025217
    headuni1 |   1.041794   .1190745     0.36   0.720      .832706    1.303383
    zapaft92 |   1.144373   .1880477     0.82   0.412     .8292678    1.579211
      yrszap |   1.017463   .0162644     1.08   0.279      .986079    1.049845
     tchload |   .9831525   .0165237    -1.01   0.312     .9512944    1.016078
  cumGOAlag1 |    .967203   .0699362    -0.46   0.645     .8394006    1.114464
   cumOTlag1 |    1.07582   .0927614     0.85   0.397     .9085437    1.273895
       avgco |   1.203184   .0611289     3.64   0.000     1.089145    1.329162
      avgcom |   .8559905   .0438663    -3.03   0.002     .7741908     .946433
 pubm_pre92h |   1.657828   .1905676     4.40   0.000     1.323408    2.076755
   pasttop2y |   2.482593   .4473648     5.05   0.000     1.743887    3.534212
  pasttop2ym |   .6499354   .1179757    -2.37   0.018     .4553663    .9276401
    fulltiun |    1.80313   .5276599     2.01   0.044     1.016097    3.199771
------------------------------------------------------------------------------

.         
. stcox male age maind_B maind_C maind_E maind_G maind_H maind_I maind_M maind_N maind_P maind_R maind_Z fac_wet fac_twt fac_lbw fac_gen f
> ac_far r_doc1 r_hdoc1 r_hl1 r_rest1 yrsrank headuni1 zapaft92 yrszap tchload cumGOAlag1 cumOTlag1 avgco avgcom pubm_pre92h pasttop2y pas
> ttop2ym fulltiun, shared(persnr) /* repeated top with frailty: table 3, dependent variable = citations */

         failure _d:  h2y == 1
   analysis time _t:  (end-origin)/365.25
             origin:  time max(zapdate,d(01oct1992))
  exit on or before:  finalobs==1
                 id:  persnr

Fitting comparison Cox model:

Estimating frailty variance:

Iteration 0:   log profile likelihood =  -3399.966  
Iteration 1:   log profile likelihood =  -3399.966  (backed up)
Iteration 2:   log profile likelihood = -3399.9652  
Iteration 3:   log profile likelihood = -3399.9652  

Fitting final Cox model:

Iteration 0:   log likelihood = -4105.0047
Iteration 1:   log likelihood = -3912.5429
Iteration 2:   log likelihood = -3618.4827
Iteration 3:   log likelihood =  -3412.684
Iteration 4:   log likelihood = -3400.0626
Iteration 5:   log likelihood = -3399.9652
Iteration 6:   log likelihood = -3399.9652
Refining estimates:
Iteration 0:   log likelihood = -3399.9652

Cox regression --
         Breslow method for ties                Number of obs      =      5284
         Gamma shared frailty                   Number of groups   =       908
Group variable: persnr

No. of subjects =          908                  Obs per group: min =         1
No. of failures =          624                                 avg =  5.819383
Time at risk    =         5284                                 max =         9

                                                Wald chi2(35)      =    305.37
Log likelihood  =   -3399.9652                  Prob > chi2        =    0.0000

------------------------------------------------------------------------------
          _t | Haz. Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        male |    6.22066   3.478021     3.27   0.001     2.079327    18.61016
         age |   .9330824   .0132988    -4.86   0.000      .907378     .959515
     maind_B |   2.457115   .9433123     2.34   0.019     1.157818    5.214477
     maind_C |   3.091652   1.155901     3.02   0.003     1.485747    6.433337
     maind_E |   2.535691   1.112486     2.12   0.034     1.073119    5.991626
     maind_G |   3.887668   2.298138     2.30   0.022     1.220432     12.3841
     maind_H |   4.525554   2.111884     3.24   0.001     1.813226    11.29514
     maind_I |    1.70228   .6391042     1.42   0.157     .8155601     3.55309
     maind_M |   2.406062   .8769641     2.41   0.016     1.177769    4.915341
     maind_N |   4.161596   2.439416     2.43   0.015     1.319193    13.12839
     maind_P |   3.324409   1.453019     2.75   0.006      1.41148    7.829863
     maind_R |   2.645032   1.318073     1.95   0.051     .9960019     7.02428
     maind_Z |   1.424451   .6049223     0.83   0.405     .6196797     3.27437
     fac_wet |   1.786717   .9889861     1.05   0.294     .6038058    5.287058
     fac_twt |   2.091231   1.184054     1.30   0.193     .6893695    6.343837
     fac_lbw |   7.141028   4.128994     3.40   0.001     2.299269    22.17847
     fac_gen |   4.739106   2.347966     3.14   0.002     1.794622    12.51468
     fac_far |   1.028767   .6956084     0.04   0.967     .2733826    3.871357
      r_doc1 |   .3204245   .0876019    -4.16   0.000      .187505    .5475684
     r_hdoc1 |   .5813765   .1156476    -2.73   0.006     .3936735    .8585759
       r_hl1 |   .6430087   .1123723    -2.53   0.012      .456521    .9056763
     r_rest1 |   .6995172   .1771467    -1.41   0.158     .4258327      1.1491
     yrsrank |   1.008247   .0197164     0.42   0.674     .9703351    1.047641
    headuni1 |   1.157556   .1510478     1.12   0.262     .8963343    1.494907
    zapaft92 |   1.014338    .224661     0.06   0.949     .6571342     1.56571
      yrszap |   1.019096   .0228497     0.84   0.399     .9752811    1.064879
     tchload |   .9795002   .0194285    -1.04   0.296     .9421518    1.018329
  cumGOAlag1 |   1.031674   .1077894     0.30   0.765     .8406375    1.266124
   cumOTlag1 |   1.224106   .1260325     1.96   0.050     1.000416    1.497814
       avgco |    1.20123   .0851899     2.59   0.010     1.045346     1.38036
      avgcom |   .8686941   .0622036    -1.97   0.049     .7549459    .9995807
 pubm_pre92h |   2.551052   .5401847     4.42   0.000     1.684526    3.863324
   pasttop2y |   2.195519    .606829     2.85   0.004     1.277233    3.774023
  pasttop2ym |   .6090311   .1692889    -1.78   0.074     .3532123     1.05013
    fulltiun |   1.756028   .4722979     2.09   0.036     1.036559    2.974877
-------------+----------------------------------------------------------------
       theta |   1.076038   .2853724
------------------------------------------------------------------------------
Likelihood-ratio test of theta=0: chibar2(01) =    35.26 Prob>=chibar2 = 0.000

Note: standard errors of hazard ratios are conditional on theta.

. 
. log close
      name:  <unnamed>
       log:  C:\Users\Stijn\Documents\_Projects_\Research\Papers\Prestatieanalyse\Persistence_in_Top\Analyse\130401_ReStat_scripts_data\su
> bmitted_to_ReStat\citanalysis.log
  log type:  text
 closed on:   8 Apr 2013, 14:01:38
------------------------------------------------------------------------------------------------------------------------------------------
